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Abstract 

The concentration of paramagnetic trace impurities in glasses can 
be determined via precise SQUID measurements of the sample’s mag¬ 
netization in a magnetic field. However the existence of quasi-ordered 
structural inhomogeneities in the disordered solid causes correlated tun¬ 
neling currents that can contribute to the magnetization, surprisingly, 
also at the higher temperatures. We show that taking into account such 
tunneling systems gives rise to a good agreement between the concen¬ 
trations extracted from SQUID magnetization and those extracted from 
low-temperature heat capacity measurements. Without suitable inclusion 
of such magnetization contribution from the tunneling currents we find 
that the concentration of paramagnetic impurities gets considerably over¬ 
estimated. This analysis represents a further positive test for the struc¬ 
tural inhomogeneity theory of the magnetic effects in the cold glasses. 


1 Introduction 

Multi-component silicate glasses are technologically important insulating struc¬ 
tural materials, usually containing trace paramagnetic impurities (typically Fe^^ 
and Fe^"*" from the fabrication process). Such dilute impurities give rise to 
nearly-ideal Langevin paramagnetism that can be exploited in low-temperature 
thermometry [T]. A by-product of the magnetization measurements (usually 
through SQUID magnetometry) is the determination of the impurity concen¬ 
tration (in the atomic ppm region). For example: commercial borosilicate 
glasses BK7 (optical glass) and Schott’s Duran (laboratory glass) have (report¬ 
edly [H m |3]) about 6 ppm and 120 ppm (or 180 ppm [1]) of diluted iron 
impurities, respectively, according to SQUID-magnetometry. 
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Knowledge of such paramagnetic impurity concentration is important also where 
the fundamental physics of disordered solids is concerned. At low temperatures 
(below 1 K, normally) the physics of a disordered solid is known to be domi¬ 
nated by low-energy excitations that go under the name of tunneling systems 
(TSs) a. These TSs are dynamical defects that give rise to quasi-universal 
physical properties that can also be exploited in low-temperature thermome¬ 
try 0. A celebrated example is the excess TS contribution to the real part e' of 
the dielectric constant at low frequency, which depends logarithmically on the 
temperature T [3I71IH]. After 40 years of research, however, the precise micro¬ 
scopic nature of the TSs is still a mystery. More recently, this dielectric constant 
e' has been discovered to be sensitive to weak magnetic fields {B ^ 10^ to 10^ 
G) [9] for some multi-component silicate glasses and possible explanations for 
this unexpected magnetic effect (observed also in other physical properties [HIS]) 
that have been proposed involve nuclear quadrupole moments nni, structural 
inhomogeneities m and also paramagnetic impurities m- Therefore, a precise 
determination of the impurity concentration is essential also in order to decide 
among different explanations for the magnetic effects. 

Research on the ultimate nature of the TSs is, moreover, receiving renewed in¬ 
terest in view of the fact that the TSs have been recognized to be the cause 
of decoherence in Josephson-j unction based quantum computing devices (the 
tunneling JJ barrier being typically amorphous) [13]. Moreover, studies of ag¬ 
ing [13] in glasses (hard and polymeric) [TS] at very low temperatures depend on 
a more complete description of the physics of the TSs. The issue of the origin 
of the magnetic effects therefore helps in improving knowledge about the na¬ 
ture of the TSs so as to minimize [T6l[T7] (or exploit [18]) their decoherence (or 
coherence) effect in SCJJ qubits. The study of the structure of real glasses at 
low temperatures provides, in turn, alternative information on the mechanism 
for the glass transition (normally investigated from the liquid state [TO] [20] ). 

In this paper we investigate mainly the issue of the determination via SQUID 
magnetization of the concentration nj of paramagnetic impurities [T], also in 
view of the fact that the theoretical analysis [U of the magnetic effect [1 
m] in the heat capacity Cp of some multi-silicate glasses produced values of 
nj systematically much lower than those quoted in the literature [Dmnig 
(and obtained from SQUID measurements). We briefly review the analysis of 
some Cp{T,B) data in the range 0.6 to 1.3 K [2], then apply our model to 
the calculation of the TS contribution to the magnetization M and analyze the 
available magnetization M (T, B) data in the range 4 to 300 K [HE]|3] with our 
formula added to Langevin’s contribution from the paramagnetic impurities. We 
find that the concentrations nj of such impurities extracted from both types of 
measurements will get to agree with each other only when the TS contributions 
(from our model, e.g.) to both Cp{T,B) and M{T,B) get to be added to 
Langevin’s known expressions for the paramagnetic impurities’ contributions. 
We apply our analysis to available data for the borosilicate Duran glass and 
for the multi-silicate glass of composition Al2 03-Ba0-Si02 (in short AlBaSiO, 
or BAS, reported concentration npe — 100 ppm [111319]) since these glasses 
have shown the most remarkable magnetic effects [3]- The concentrations nj 


2 


that we find for the Fe impurities are, however, about 60 to 80% lower than 
those quoted in the literature for these glasses, so the correct description of 
the magnetic-sensitive TSs becomes important also for the applications of low- 
temperature physics. 

This paper is organized as follows: in Section 2 we present the foundations of 
the extended tunneling model that we use, then in Section 3 we apply the model 
to the explanation of the magnetic effect in the Cp{T, B) data. In Section 4 we 
derive the contribution from the TSs in our model to the magnetization M{T, B) 
and analyze the available data with our formula. The values of nj extracted 
from both quantities Cp and M are compared in Section 5, which also contains 
our conclusions. In the Appendix we present some preliminary results for the 
analysis of some measurements in BK7 glass. 


2 The Extended Tunneling Model 


The modern justification for the tunneling model is based on the belief that 
glasses at sufficiently low temperatures are characterized by a potential energy 
landscape (PEL), that can be investigated only by means of classical computer 
simulations of atomic configurations. Out of the many local minima of the PEL 
some local potentials, normally thought to be double-welled (DWP), give rise 
to TSs, normally replaced by two-level systems (2LSs) that at low temperatures 
are characterized by the tunneling Hamiltonian [B]: 


'^2LS 


1 / A Ao \ 

2 Ao -A ; • 


( 1 ) 


Here the parameters A (the energy asymmetry) and Aq (twice the tunneling 
parameter) are typically characterized by a probability distribution that views 
A and In(Ao) (the latter linked to the DWP energy barrier) broadly (in fact 
uniformly) distributed throughout the disordered solid [^: 

'P2Ls(A, Ao) = — (2) 

where some cutoffs are introduced when needed and where P is a material- 
dependent parameter, like the cutoffs. In this justification of the TSs clearly the 
tunneling “particle” cannot be typically a real atom/ion of the glass, but rather 
an effective, fictitious particle: jumps between contiguous low-lying minima of 
the PEL correspond to the rearrangements of several atoms/ions, if not large 
parts of the entire atomic configuration. 

Much progress has been done in understanding the physics of glasses at low 
temperatures with the help of the above simple model (the standard tunnel¬ 
ing model, STM) [3], which has however important limitations. The recently 
discovered magnetic effects in non-magnetic glasses mil II [To] , for example, 
cannot be explained without a suitable extension of the STM. Our own ex¬ 
tension brings three realistic considerations in the modeling of the structure of 
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real glasses. The first is that glasses can be no longer considered fully homoge¬ 
neously disordered solids at the intermediate atomic scales, for there is mount¬ 
ing experimental evidence that the structure is spatially inhomogeneous with 
“better ordered” regions being inter-dispersed in an otherwise featureless ho¬ 
mogeneously disordered matrix. One way to look at these regions of enhanced 
atomic ordering is that they are the thermal-history continuation to temper¬ 
atures below the glass transition temperature Tg [22] of the slower-particles 
regions present, within the sea of faster-particle regions, in the dynamical- 
heterogeneity picture [53] of the supercooled liquid phase, between Tg and the 
melting temperature Tm > Tg. These slower-particle regions are in fact also 
better ordered, as expected. We name these regions in the glassy phase regions 
of enhanced regularity (RERs), but some other names have been proposed in 
the literature: cybotactic groupings from a critical analysis of the X-ray and 
neutron scattering data in amorphous solids |24] . para-crystals from combined 
electron-diffraction and fluctuation electron microscopy of a-Si films |25j . Simi¬ 
lar conclusions about partial devitrification have been reported for the metallic 
glass Zr 5 oCu 45 Al 5 using combined Monte Carlo simulation and fluctuation elec¬ 
tron microscopy |26j . Besides, the amorphous solids of general composition 
(Mg 0 )a;(Al 203 )y(Si 02 )i-x-y are termed ceramic glasses and are even known 
to contain embedded micro-crystals m- The evidence from X-ray analysis for 
crystalline-like ordering in quenched network glasses is in fact most compelling 
in the case of multi-component materials [33] . For these glasses the distribution 
([2]) of the 2LS parameters should be partially abandoned in favour of a different 
distribution, for a subsystem of TSs nesting within the RERs, having a form 
favouring the near-symmetry of the local wells of the DWP, e.g.: 

= (3) 

where the parameters A*, Aq and P* now refer to that subsystem. One should 
then use at the very least m a collection of the two types of TSs described above 
when dealing with multi-component glasses (for which the strongest magnetic 
effects have been observed). The second consideration is that the computer¬ 
generated PEL typically contains more complicated, local multi-welled poten¬ 
tials (MWPs) as well as DWPs (indeed, a representation of the PEL solely by 
means of DWPs and 2LSs seems rather oversimplified, though of obvious prac¬ 
tical theoretical advantage [IH]). In Fig. |T] the one-dimensional map of a molec¬ 
ular dynamics-generated PEL for a system of 32 particles interacting through 
a binary mixture Lennard-Jones potential in the glassy phase has been repro¬ 
duced [29]. There are clearly DWPs, but also MWPs as we have highlighted. 
The tunneling Hamiltonian of a MWPs is easily written down as m- 


El 

Do 

Do \ 


Do 

E 2 

Do 

(4) 

Do 

Do 

Eo / 



where Ei, E 2 , E 3 are the energy asymmetries between the wells (we have chosen 
the simplest case, having three wells) and Dq is the most relevant tunneling 
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Figure 1: The PEL for a BLJ mixture of 32 particles, mapped in ID according 
to the procedure by Heuer. Highlighted in black are the 2LSs, then in blue the 
3LSs, 4LSs, ... ensuing from multi-welled local potentials. Adapted from pU] 


amplitude (through saddles of the PEL, in fact). This 3LS Hamiltonian has the 
advantage of readily allowing for the inclusion of a magnetic field B > 0, when 
coupling orbitally with a tunneling “particle” having charge q {q being some 
multiple of the electron’s charge —e) m-- 

/ El L>oe*^/3 Doe-*^/3 \ 

nsLsiB) = E 2 Doe^^/^ (5) 

\ Doe-*‘^/3 E 3 ) 

where i^/3 is the Peierls phase for the tunneling particle through a saddle in the 
field, and (p is the Aharonov-Bohm phase for a tunneling loop and is given by 
the usual formula: 

$ h 

iP = 2tt—, ^0 =-tt (6) 

^0 \q\ 

$0 being the appropriate flux quantum {h is Planck’s constant) and $ = B • S 
the magnetic flux threading the area S formed by the tunneling paths of the 
particle in this simple model. The energy asymmetries Ei, E 2 , E 3 typically 
enter through their combination D = i/Ef + E 2 + E^. One can easily convince 
oneself that if such a MWP is used with the standard parameter distribution, 
Eq. ([2) with D,Dq replacing A,Ao, for the description of the TS, one would 
then obtain essentially the same physics as for the STM 2LS-description. In 
other words, there is no need to complicate the minimal 2LS-description in 
order to study glasses at low temperatures, unless structural inhomogeneities of 
the RER-type and a magnetic field are present. Without the RERs, hence no 
distribution of the type ([3]), the interference from separate tunneling paths is 
only likely to give rise to a very weak Aharonov-Bohm effect. Hence, it will be 
those TSs nesting within the RERs that will give rise to an enhanced A-B effect 
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and these TSs can be minimally described - for example - through Hamiltonian 
([5]) and with distribution ([3]) thus modified to favour near-degeneracy m- 


r;i^siEuE2,E3;Do) = 


p* 


DoiEf 


El 


EIY 


(7) 


We remark that the incipient “crystallinity” of the RERs calls for near-degeneracy 
in El, E 2 , E^ simultaneously and not in a single one of them, hence the corre¬ 
lated form of ©• Other descriptions, with four-welled potentials or modified 
three-dimensional DWPs are possible for the TSs nested in the RERs and lead 
to the same physics as from Eqs. and © above [30] (which describe what 
we call the anomalous tunneling systems, or ATSs, nesting within the RERs). 
The final and most important consideration is that the TSs appear to be rather 
diluted defects in the glass (indeed their concentration is of the order of mag¬ 
nitude of that for trace paramagnetic impurities, as we shall see), hence the 
tunneling “particles” are embedded in a medium otherwise characterized only 
by simple acoustic-phonon degrees of freedom. This embedding, however, means 
that the rest of the material takes a part in the making of the tunneling poten¬ 
tial for the TS’s “particle”, which itself is not moving quantum-mechanically in 
a vacuum. Sussmann m has shown that this leads to local trapping poten¬ 
tials that (for the case of triangular and tetrahedral perfect symmetry) must be 
characterized by a degenerate ground state. This means that, as a consequence 
of this TS embedding, our minimal model ([5]) must be chosen with a positive 
tunneling parameter m-- 

Do>0 ( 8 ) 

where of course perfect degeneracy is always removed by weak disorder in the 
asymmetries. The intrinsic near-degeneracy of 0 implies that this model should 
be used in its D/Dq <§; 1 limit, which in turn reduces the ATSs to effective 
magnetic-field dependent 2LSs and greatly simplifies the analysis together with 
the limit </?—>■ 0 which we always take for relatively weak magnetic fields. Our 
extended tunneling model (ETM) consists then in a collection of independent, 
non-interacting 2LSs described by the STM and 3LSs described by Eqs. ([5]) and 
© above in the said D/Dq <C 1 and ip ^ 0 limits, the 3LSs nested within 
the RERs and the magnetic-field insensitive 2LSs distributed in the remaining 
homogeneously-disordered matrix. Our ETM has been able to explain the mag¬ 
netic effects in the heat capacity m, in the real |32] and imaginary |33j parts 
of the dielectric constant and in the polarization echo amplitude |33j measure¬ 
ments reported to date for various glasses at low temperatures, as well as the 
composition-dependent anomalies [8l|3l]. The new physics is provided by the 
magnetic-field dependent TS density of states (DOS) which acquires a term due 
to the near-degenerate MWPs m that gets added up to the (nearly) constant 
DOS from the STM 2LSs (having density n 2 Ls)- 


gtot{E, B) = n2LsP + nATS-prfATs{E, B)9{E — Ec) (9) 

h/ 

where uats is the ATSs’ concentration, Jats is a magnetic-field dependent 
dimensionless function, already described in previous papers m, and Ec is a 
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Figure 2: The magnetic-sensitive part of the density of states (DOS) as a func¬ 
tion of the A-B phase ip (proportional to the magnetic field B) and different 
energies {uatsP* has been set to 1). The shape of this part of the DOS (coming 
from the MWPs with a parameter distribution 0 favouring near-degeneracy) 
is the ultimate source of all the magnetic effects. The cusp is an artifact of the 
effective 2LS approximation but also of the existence of upper and lower 
bounds for Dq owing to the nature of the RER atomic structure. 


material and R-dependent cutoff. The 1/E dependence is a consequence of the 
chosen tunneling parameter distribution, Eq. o (or, to the same effect, Eq. ([3])), 
and gives rise to a peak in gtot near Ec that is rapidly eroded away as soon as 
a weak magnetic field is switched on. The form and evolution of the magnetic 
part of the DOS is shown in Fig. [5] for some typical parameters, as a function 
of B for different values of E. This behaviour of the DOS with B is, essentially, 
the underlying mechanism for all of the experimentally observed magnetic field 
effects in the cold glasses within this model: the measured physical properties 
are convolutions of this DOS (with appropriate B-independent functions) and 
in turn reproduce its shape as functions of B. As an example, the total TS heat 
capacity is given by 

pOO 

CpTsiT,B)= dE gtot{E,B)CpoiE,T) 

Jo 

where 

Cpo{E,T) = kB cosh 2 

is the heat capacity contribution from a single TS having energy gap E. 
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We end this Section by mentioning that in model simulation studies of the 2LSs 
generated by defects in perfect crystals, Churkin, Barash, and Schechter [35] 
found non-uniform behaviour for the g{E) DOS very similar to what we advocate 
in Eq. (jH]) for the ATS part (see Ref. [TT], Fig. 6 for B=0). The (near-) 
crystallinity is thus the common theme between the RERs in glasses and the 
twin-caged defects in crystals. 

A much deeper justification for our ATS model and for the nature of the TSs 
in general (of two types only: the 2LSs and the three- or four-fold ATSs) will 
be presented elsewhere, for glasses (papers in preparation). 


3 Heat Capacity 


3.1 Summary of previous work [ 111] 


In this Section we re-analyze the available data [5] for the magnetic effect in the 
heat capacity of two multi-component glasses, commercial borosilicate Duran 
and barium-allumo-silicate (BAS) glass, in order to better estimate the concen¬ 
tration of trace Fe-impurities in this way. The presence of a magnetic effect in 
the Pyrex glass was reported long ago by Stephens m and attributed solely to 
paramagnetic iron impurities even though the maximum effect was for R ~ 0. 
A systematic experimental study of Cp{T,B) around and below 1 K in some 
multi-silicate glasses was carried out by Siebert [2] and those data have been 
used, upon permission, by one of us m as the very first test of the above (Sec¬ 
tion 2) ETM [11]. That earlier analysis best-fitted the Cp data by Siebert with 
the sum of Einstein’s ^phT^ phonon term plus the 2LS ^ 2 LsT non-magnetic 
contributions, as well as with Langevin’s paramagnetic and the ATS contribu¬ 
tions (see below). The analysis came up with concentrations nj ~48 ppm and, 
respectively, fij ~ 20 ppm instead of the quoted [5] 126 ppm (or 180 ppm in a 
different study DP) and 102 ppm for Duran and for BAS glass, respectively. 

In order to better understand this large discrepancy we begin with by re¬ 
analysing Siebert’s data for Cp{T,B), after subtraction of the data taken at 
the same temperatures for the same glass, but in the presence of the strongest 
applied magnetic field (8 T). In this way only the magnetic-field dependent con¬ 
tributions should remain in the data for Cp{T,B) = Cp{T,B) — Cp{T,oo). It 
should be remarked that use of these data for Cp is marred by the question 
of the many relaxation times in a glass, since the data are obtained indirectly 
(like in all Cp measurements in glasses at low T), from a heat-pulse experiment 
where the change with time of the sample’s temperature is fitted with a single 
relaxation time involving Cp. The data thus obtained for Cp, however, strongly 
resemble those obtained earlier on by Stephens (at B=0 and 3.3 T only) [21] and 
there is theoretical evidence that the width of the relaxation times’ distribution 
narrows considerably in a magnetic field for a multi-component glass [36] . Still, 
we use the data of [2] tentatively. 

Fig. [3] presents the data for the heat capacity after subtraction of the mag¬ 


netic field independent data at B=8 T for the BAS glass (Fig. 3(a)) and for 
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Duran (Fig. |3(b)[ ). The parameters to be determined are the cutoff Dmin and 
combinations of cutoffs, charge and area DominQS and Domax^S [a, as well as: 

np^ 2 + Fe^'*" impurity concentration 
np^3+ Fe^"*" impurity concentration 
nATS ATS concentration (always multiplied by P*) 


The data of Fig. |3(a)| and Fig. |3(b)| (we have restricted the best fit to the three 
temperatures having the most data points around the peak of Cp{B)) have been 
best-fitted by using the following magnetic-dependent contributions: 

a. the known Langevin contribution of the paramagnetic Fe impurities (Fe^'*" 
and Fe^’*') having concentration nj: 


( 10 ) 


where z = and where g is Lande’s factor for the paramagnetic ion 

in that medium, /tb is Bohr’s magneton and J the ion’s total angular mo¬ 
mentum (in units h=l)\ kp is Boltzmann’s constant. We have assumed 
the same values the parameters g and J take for Fe^+ and Fe^’*' in crys¬ 
talline Si02: J=2 with g=2 and, respectively, J = 5/2 with g=2 (we 
have adopted, in other words, complete quenching of the orbital angular 
momentum m, consistent with other Authors’ analyses mm)- 

b. the averaged contribution of the ATSs written in terms of a sum of 
individual contributions from each ATS of lowest energy gap E 


Cats{T, (f) = 



TT P* UATS 
4 kpT^ 

E 


cosh^( 


E 

2kBT 


) 


-b 



dE 


E 

(ikpr) 


In 


In 




DminDormn'f^ 

/ Domax \ ^ ^ ~ ^Omin^P 
V n„ . 1 


- Domin E'^ ^OmaxP^- 


( 11 ) 


or, re-written in a dimensionless form as 


Cats{T, ip) = Co{T, p) + 27rP* nATsfesj [l^( ^cl) -^{^c2)] 

+ i \T{Xc 1 i^min) H” 7 ^Omin p) -X{Xc2,Xo max 

( 12 ) 


where: 


• Ed = VE'Lu + and E ^2 = 


• Xd,2 


Ec1,2 . 

2feBT’ 


2kBT^ 


etc.; 
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BAS glass 



X 10"^ 


(a) 

Duran 



(b) 


Figure 3: The heat capacity data Cp = Cp — Cp{8 T) as a function of the 
magnetic field B for different temperatures: a) for the BAS and b) Duran 
glasses. Data from [5] (also reproduced, upon permission, in [TT]1. 
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• I{x) = X tanh a: — In cosh x] 

• ^( 2 ^- «) = XT ln(y2 _ a^). 

and the following expression: 


Co(T, = 27rP*nATskB In ( 11^(2) - /(x^z).} (13) 

\ ^Omin J 

The angular average over the ATS orientations is performed by replacing 
(/5 —(in other words averaging cos^ 9, 9 being the orientation of S 
with respect to B). 

3.2 Extracted parameters from best-fitting the heat ca¬ 
pacity data 

3.2.1 BAS glass 

The concentrations of the ATSs and Fe-impurities extracted from the best fit of 
the heat capacity as a function of B, for the BAS glass, are reported in Tabled! 
having fixed the concentrations, it was possible to extract the other parameters 
for the BAS glass (Table dj. The best fit of the chosen data is reported in 
Fig-M 


BAS glass 

Concentration [g 

Concentration [ppm] 

TIpq2+ 

1.06x10^'' 

14.23 

nf'Q3+ 

5.00x10^® 

6.69 

P*nATS 

5.19x101® 

- 


Table 1: Extracted parameters (from the heat capacity data) for the concentra¬ 
tions of ATSs and Fe-impurities for the BAS glass. 


Temperature [K] 

Dmin [K] 

Dormn\i\S [KA^j 

Doraa.\i\S [KA^] 

0.60 

0.49 

4.77x10-1 

3.09x10® 

0.90 

0.53 

5.07x10-1 

2.90x10® 

1.36 

0.55 

5.95x10-1 

2.61x10® 


Table 2: Extracted tunneling parameters (from the Cp data) for the BAS glass. 


3.2.2 Duran 

The concentrations of the ATSs and Fe-impurities extracted from the best fit 
of the heat capacity as a function of B, for Duran, are reported in Table |3l 
having fixed the concentrations, it was possible to extract the other parameters 
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BAS glass 



(a) 


Duran 



(b) 


Figure 4: The heat capacity best fit for the a) BAS and b) Duran glasses. 
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for Duran (Table S]). The fit of the chosen data is reported in Fig. |4(b)[ The 
first comment we make is that these good fits, with a smaller set of fitting 
parameters, repropose concentrations and tunneling parameters very much in 
agreement with those previously obtained by one of us m- The problem with 
the concentrations of the Fe-impurities reported in the literature is that they 
do not allow for a good fit of the Cp{B) = Cp{B) — Cp{oo) data in the small 
field region when the Langevin contribution alone is employed (Eq. (jlO[) without 
Eq. (fT^ 'l. See for instance Section 5, Eig. [T] The Langevin contribution drops 
to zero below the peak, whilst both Siebert’s and Stephens’ data definitely point 
to a non-zero value of C'p(O) = C'p(O) — Cp{oo) at B — 0 for any T > 0. This 
non-zero difference is well accounted for by our ETM and it comes from the 
ATS contribution to the DOS in Eq. (l9|). 


Duran 

Concentration [g 

Concentration [ppm] 

TIpq2+ 

3.21x10^^ 

33.01 

TIpq3+ 

2 .11x10^^ 

21.63 

P*nATS 

8 .88x10^6 

- 


Table 3: Extracted parameters (from the heat capacity data) for the concentra¬ 
tion of ATSs and Fe-impurities for Duran. 


Temperature [K] 

Dmin [K] 

Dormn\i\S [KA^j 

Doma.\i\S [KA^j 

1.11 

0.34 

4.99x10-^ 

2 .68x10^ 

1.23 

0.32 

5.30x104 

2.50x10® 

1.36 

0.32 

5.54x104 

2.46x10® 


Table 4: Extracted tunneling parameters (from the Cp data) for Duran. 


The results of our Cp analysis definitely indicate that the concentration of 
paramagnetic impurities in the multi-silicate glasses is much lower than previ¬ 
ously thought and extracted from SQUID-magnetometry measurements of the 
magnetization M{T,B) at moderate to strong field values and as a function 
of T. Therefore we now turn our attention to a re-analysis of the SQUID- 
magnetometry data. 

4 Magnetization 

4.1 Theory 

The second comment we make, having obtained a qualitatively (and also quan¬ 
titatively, barring the multiple relaxation-times issue) good fit to the Cp{B) 
data with the ATS contribution added to the Langevin’s, is that the ATSs now 
appear to carry considerably high magnetic moments (J-ats- Estimating from 
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the definition (T = 0) (J-ats = ~'§b{~ 1^)^ where E = ^JD'^ + D^ip^ is the 
ATS lower energy gap, we get for not too small fields B (hats vanishes linearly 
with B when B —5> 0, but saturates at high enough B): 

^^^TS-^SDQ = ^(\^SDo). (14) 

$0 00 ^e' / 

Thus the very same combination ^SDq of parameters appears, whilst (j^o = h/e 
is the electronic magnetic flux quantum. Using the values extracted from the 
Cp best fit (e.g. TablelU we deduce from Eq. (ITTl) that (for Duran) fiATS ranges 
from about 3.8^b to 27.1^b- This fact alone indicates that a large group of 
correlated charged atomic particles is involved in each single ATS and that an 
important ATS contribution to the sample’s magnetization is to be expected 
(Fe^'*' and Fe^"'' have magnetic moment p,j = A^b and, respectively, 5fiB)- 
The magnetization M of a sample containing dilute paramagnetic impurities as 
well as dilute magnetic-field sensitive ATSs is, like Cp, also given by the sum of 
two different contributions: 


a. Langevin’s well-known paramagnetic impurities’ contribution (Fe^'*" and 
Fe^“'", with nj concentration of one species having spin J), given by the 
standard expression [35] 


Mj = njgfj.BJBj{z), 


z = 


ggBBJ \ 
kBT ) 


where the Brillouin function Bj is defined by: 


Bj{z) = 


2J- 


2J 




2J 




(15) 


(16) 


and its low-field susceptibility is the known Curie law: 

^ ^ njg'^g.BJjJ + 1 ) 

B ~ 3kBT 


(17) 


b. the ATS tunneling currents’ contribution, given by the following novel 
expression as the sum of contributions from ATSs of lowest gap E: 

C^c2 


Mats = P*nATS-fz 

ID 


' Eel 
nCO 

fEe2 


\2kBT) V DCn ) 


min 
2 


dE tanh ( —In (^ 
\2kBTj 


(18) 


and which can be also re-expressed (like in the case of Cats) using y = 
2 j^rp in the following form: 

Mats = 27r P*nATskBT— 

ID 


-L 


Xcl 

CO 


dy tanh y In 
dy tanh y In 






(19) 
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with, as before: 


• Ed = V+ Dl^in^P^ and Ec 2 = 


• Xd,2 


Eel,2 
2 /cbT’ 


^min 


Dmin 

2fcBT’ 


etc.; 


We present this expression here for the first time, also motivated by the fact that 
we expect a contribution to the measured magnetization M from the ATSs that 
is comparable to, or even greater than, Langevin’s paramagnetism of the diluted 
Fe impurities. The above expression follows from a straightforward application 
of standard quantum statistical mechanics, with 


Mats = nATs{ 


dHsLS V 

dB 


nATS being the ATSs’ concentration (a parameter always lumped together with 
P*) and with 'H^ls given by Eq. ([5]). The angular brackets (• • •) denote quan¬ 
tum, statistical and disorder averaging. 

The above formula for Mats is in fact correct for weak magnetic fields. For 
higher fields a correction has to be introduced, owing to the fact that an im¬ 
proved analytic expression [39] for the lowest ATS energy gap must be used: 


e = ^d^ + d1^^i-E^^). 

A full derivation and the study of the B- and T-dependence of Mats will be 
presented elsewhere. 


4.2 Extracted parameters for the magnetization data 

4.2.1 BAS glass 

The magnetization data [5] were best-fitted with Eq. (lisp (for the Ee^"*" and Fe^+ 
contributions) as well as with Eq. (fT^ (for the ATSs’), using the parameters 
from the Cp-fits as input. The best fit for the BAS glass is reported in Fig. [5] 
and the extracted parameters in Table jSj 

4.2.2 Duran 

The best fit for Duran is reported in Fig. [6] and the extracted parameters in 
Table El 


4.3 Concentration conversion 


We now explain how we convert the Fe-concentrations thus obtained to atomic 
ppm concentrations (ppma). nj, the mass density of a Fe species with spin J 
in the sample, is given by: 


Nj Nj Na 
M Nat E. 


( 20 ) 
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BAS glass (B = 3 T) 



Figure 5: The best fit of the magnetization data [2] for the BAS glass, using 
Eq. (fTK]) (for the Fe^+ and Fe^+ impurities) and Eq. ([T^ (for the ATSs). 


Parameter 

BAS glass 

«Fe2+ [g"-^] 

1.08x10^'' 

«Fe3+ [g"^] 

S.OlxlQi® 

P*nATS [g“^] 

5.74x101® 

Dmin [K] 

8.01x10-2 

[Kl2] 

1.31x10® 

Doraax\l\S [KA^] 

2.44x10® 

vert.offset [Am^g“^] 

-1.04x10"® 


Table 5: Extracted parameters (from the magnetization data of [5]) for the 
concentration of ATS and Ee-impurities of the BAS glass. The vertical offset 
represents the residual Larmor diamagnetic contribution. 


where 




is the sample’s mass, and where: 


( 21 ) 


16 











Duran (_B = 5 T) 



Figure 6: The best fit of the magnetization data [5] of Duran, using Eq. m 
(for the Fe^+ and Fe^'*' impurities) and Eq. (ITO)) (for the ATSs). 


Parameter 

Duran 

nFe^+ [g“-^] 

«Fe3+ [g"^] 
P*nATS [g“^] 
Pmin [K] 
Dt,rmn\l\S [KA^] 

[KAA 

vert.offset [Am^g“^] 

3.07x10^^ 
2.13x101^ 
8.68x101® 
5.35x10-2 
2.00 xlO® 
2.81x10® 
-1.97x10"® 


Table 6: Extracted parameters (from the magnetization data of [5]) for the 
concentration of ATS and Fe impurities of Duran. 


Nj number of Fe-ions in the sample with spin J 
molar fraction of the i-th species 
Ai molar mass of the i-th species 
Nat total number of atoms in the sample 
Na_ Avogadro’s number (6.022x10^^ mol“^) 
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Oxide Element 

Molar Mass (A^) 

BAS glass 

6.% Duran 

Si02 

60.084 

72.7 

83.4 

B2O3 

69.620 

0.72 

11.6 

AI2O3 

101.961 

00 

bo 

1.14 

Na20 

61.979 

0.28 

3.4 

K2O 

94.196 

0.064 

0.41 

BaO 

153.326 

17.0 

0.005 

Li20 

29.881 

0.014 

0.004 

PbO 

223.199 

0.48 

<0.01 


Table 7: Molar mass Ai and percentage fraction % of the various oxides mak¬ 
ing up the BAS (third column) and Duran (fourth column) glasses, as reported 

in [3113]. 


for the Fe^+ (■^=2) and Fe^’*' (J=5/2) impurities. Tabled shows parameters 
related to the chemical and molar composition |2] of the two multi-component 
silicate glasses. Therefore, using the parameters reported in Tabled one has: 




80.530^ for BAS glass 
61.873-2-7 for Duran 


We interpret Nj/Nat = nj as the atomic concentration of the spin-J Fe species 
(to be multiplied by 10® to obtain the ppm) and thus we have the conversion 
formula: 


nj = nj 


Na 


( 22 ) 


5 Conclusions: Extracted Concentrations for Iron 
Impurities and ATSs 

5.1 BAS glass 

The nominal concentration of Fe®’*' for the BAS glass is (using Eq. (12(111 or (l22|) 1: 
fip^3+ = 102 ppm 

10-® • 102 ■ 6.022 X 10^3 mor^ „ .7.17-1 (23) 

"*•4;. =-sassoi-= 

which is inadequate (as in Duran’s case below) to explain the behaviour of the 
heat capacity as a function of B, here presented (Fig. |4|) and as a function of T 
(studied in |11]1. Table |8| summarizes the concentrations found from our best 
fits of heat capacity and magnetization data, for the BAS glass. The parameters 
found in m were P*n74TS=6.39x 10^® g ^ and nFe=20.44 ppm where this latter 
was for the Fe^"*- concentration only. The present study confirms that most of 
the Fe-impurities in these two glasses are of the Fe^+ ^ypg |Ti] . 
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BAS glass 

Heat Capacity fit 

Tlp^2+ 

1.06x10^'^ g ^ = 14.23 ppm 

npf.3+ 

5.00x10^® g“^ = 6.69 ppm 

P*nATS 

5.19x101® g-i 

Magnetization fit 

npf.2+ 

l.OSxlQii g-i = 14.38 ppm 

npf,3+ 

5.01x101® g“i = 6.70 ppm 

P*nATS 

5.74x101® g-i 


Table 8: Comparison between the concentrations extracted from the two differ¬ 
ent best-htted experimental data sets for the BAS glass. 


5.2 Duran 

The nominal concentration of Fe^"*" for Duran is (using Eq. or ([H])): 


^TLOm 


126 ppm 

10"® • 126 ■ 6.022 X 10^3 


61 . 873 ^ 

mol 


1.23- lO^V^ 


(24) 


which again is inadequate to explain the behaviour of the heat capacity as a 
function of B (see Fig. El). Table[9]summarizes the concentrations found from our 
best fits of heat capacity and magnetization data, for Duran. The parameters 
found in m were P*nA7'5=6.92x 10^® g ^ and npe=47.62 ppm where this 
latter was the Fe^’*' concentration only. 

Fig. El presents the behaviour of the different contributions to the heat capacity 
as a function of B for Duran; Cpe'^+ and Cpe3+ are given by Eq. m, respectively 
with the Fe^+ and Fe^+ parameters, Cparam is the sum of these latter two. 
Cats is given by Eq. (ini) and the green line represents the result of the best 
fit. The dashed line corresponds to the Cp{B) one would get from the nominal 
concentration fipe of 126 ppm [5] as extracted from the SQUID magnetization 
measurements fitted with the Langevin contribution only (no ATS contribution). 
Likewise Fig. [8] presents the behaviour of the different contributions lEa. dTKI) 
and Ea. (|T^ l to the magnetization as a function of B, also for Duran. It can be 
seen that the ATS contribution is in both cases dominant, also (in the case of 
the magnetization) at the higher temperatures. 


5.3 Conclusions and outlook 

In conclusion, by allowing a contribution to the magnetization from the ATS 
tunneling currents as evaluated from our ETM we have been able to show that 
the concentrations of Fe impurities (of type Fe^’*' as well as Fe^+) are in much 
better agreement when extracted from SQUID magnetization data or from heat 
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Duran (T = 1.23 K) 



Figure 7: The curves represent the different terms that contribute to the heat 
capacity of Duran in our best fit of the data from [5] . The dashed curve is for 
Langevin’s contribution only, but with the nominal concentration of npg3+=\2Q 
ppm (no ATS). 


Duran (B = 5 T) 



Figure 8: The curves represent the different contributions to the magnetization 
of Duran in our best fit of the data from [2] . 
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Duran 

Heat Capacity fit 

Tlp^2+ 

3.21x10^^ g ^ = 33.01 ppm 

nFe3+ 

2.11x10^^ g“^ = 21.63 ppm 

P*nATS 

8 .88x101® g-i 

Magnetization fit 

nFe^+ 

3.07x1011 g-i = 31.58 ppm 

nFe^+ 

2.13x1011' g“i = 21.86 ppm 

P*nATS 

8 .68x101® g-i 


Table 9: Comparison between the concentrations extracted from the two differ¬ 
ent best-htted experimental data sets for Duran. 


capacity data with the samples in a magnetic field. We have aimed at a semi- 
quantitative agreement, but inclusion of a possible contribution to Cp from 
Kramers-doublets’ splitting for the J = 5 j2 Fe^+ species [2] (which we con¬ 
firm, however, to be the minority one m) could possibly further improve the 
quantitative agreement. 

We reach the surprising conclusion that the ATSs and the ETM employed keep 
giving a good description of the glass magnetization till the highest temperatures 
in the SQUID-magnetometry measurements, about 300 K (the glass-transition 
temperature Tg being, however, 1123 K for BAS glass and 803 K for Duran, 
respectively ID). This new finding is interpreted with the consideration that 
no phonon-TS interactions are involved in the measurement of the magneti¬ 
zation, contrary to the case of AC dielectric-constants and polarization-echo 
measurements, or to the case of the acoustic measurements, or that of the heat 
capacity, where phonons do contribute and in a complicated way well above 1 
K [4l[^. The SQUID magnetization measurement in a magnetic field is thus an 
ideal arena where to study the TSs in general and to test our ETM, proposed 
for the explanation of the magnetic effects in the multi-component glasses as 
manifestations of the inhomogeneous structure of the glassy state. Therefore 
a systematic experimental study of SQUID magnetization as a function of T 
and B from multi-component glasses with very low paramagnetic impurity con¬ 
centration (e.g. BK7) would be most welcome to test our model. Moreover, 
also the tunneling parameters and the concentration of the ATSs here obtained 
turn out to be similar when extracted from the Cp- and from the M-data, with 
the tunneling parameters Dmin, Domin\^\S and Dornax\^\S extracted from the 
M-data being as anomalously large as from best-fits of all data in the other 
experiments [n [31 |33l |34]. These large values have been interpreted by one 
of us [8] as deriving from the correlated tunneling of a large, but not yet meso¬ 
scopic or macroscopic, number N ~ 200 to 600 of atomic-scale TSs, the ATS 
being only a fictitious tunneling particle involving in fact the correlated rear¬ 
rangement of a large group of (charged) atoms. There appears to be a weak 
temperature dependence N{T) of this number of correlated atomic tunnelers, 
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since the parameters quoted above change slightly (or even significantly, yet 
remaining large) from experiment to experiment carried out at different tem¬ 
perature ranges. About this revealing T-dependence will be expanded in future 
publications. 

We remark that the concentration nj of paramagnetic impurities turns out to 
be significantly lower, up to 80% less, than the concentrations reported in the 
literature for the analysed glasses and as extracted from SQUID-measurements 
of the magnetization [11 mg. Our point of view is that without the inclusion 
of the contribution from the ATSs the extracted SQUID-measurement concen¬ 
tration of paramagnetic impurities will be considerably overestimated. This 
happened already in the case of Stephens’ data jH], which have led to the 
Langevin-only estimate m npf, ~ 50 ppm in Pyrex from the Cp measure¬ 
ment data in a magnetic field, when in fact the mass-spectrometry analysis had 
given m nFe=12 ppm. The fact that we have established, that Langevin- 
only fitted SQUID-magnetization measurements considerably overestimate the 
concentration of paramagnetic impurities in a glassy matrix, would appear to 
cast serious doubts about the trace Fe-impurities and associated paramagnetic- 
TSs as possible sources of the magnetic effects in the cold glasses [ig. In¬ 
deed, Fe^+ would enter substitutionally to 81"^+ only in a crystal (e.g. quartz), 
whilst in a multi-silicate glass the overwhelming majority of Fe-impurities would 
enter as network-modifiers of the Si04 glassy matrix [ID]. This means that 
the amplitude of the [Fe04]“ paramagnetic-TS contribution should be reduced 
considerably more than the 80% we claim from the overestimate of the Fe- 
concentration from Langevin-fitted SQUID-measurements. We note in passing 
that the paramagnetic-TS explanation of the magnetic effect in Cp requires 
concentrations npe already some 40% greater than the nominal, Langevin-only 
SQUID-extracted values m- The paramagnetic-TS approach may nevertheless 
retain some validity in the case of the heavily Fe- and Cr-doped multi-silicate 
glasses [41]. 
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Appendix 

Here we present our preliminary study of the SQUID magnetization data (also 
available from m) for the borosilicate glass BK7, for which however no substan¬ 
tial magnetic effect in the heat capacity Cp has been reported [2]- This glass 
has a nominal Fe-impurity concentration of npg3+ =6 ppm niiig, yet our best 
fit in Fig. lAll with both Langevin (Eq. (HU) and ATS (Eq. (fTOl) ') contributions 
produces the concentrations and parameters given in Table lAT] The best fit was 
carried out with knowledge of ATS parameters from our own theory [33] for the 
magnetic effect in the polarization-echo experiments at mK temperatures [3]. 
We conclude that our main contention is once more confirmed, in that the con¬ 
centration of Fe in BK7 we extract in this way is only about 1.1 ppm and the 
bulk of the SQUID magnetization is due to the ATSs. Table lATl reports our very 
hrst estimate of uatsP* for BK7. Assuming P* to be of order 1 and about the 
same for all glasses, we conclude that the concentration hats of the ATSs nest¬ 
ing in the RERs is very similar for all of the multi-silicate glasses by us studied 
for their remarkable magnetic effects. From the present SQUID-magnetization 
best fits we have obtained 5.74x10^® g“^ (BAS glass), 8.68x10^® g“^ (Duran) 
and 1.40x10^® g“^ (BK7). The almost negligible magnetic effect in Cp for BK7 
is due, in our approach, to the low values of the cutoffs D^rnin and D^max for 
this system (these parameters appearing in the prefactor and in the integrals’ 
bounds determining the ATS contribution to Cp i)- 


Parameter 

BK7 

nFe2+ [g“-^] 
nFe3+ [g“^] 
P*nATS [g“^] 
Pmin [U] 
Dormnl^JS [KA^] 
Domaxl^JS [KA2] 
vert.offset [Am^g“^] 

6.69x10^® = 0,71 ppm 

3.43x10^® = 0.36 ppm 
1.40x101® 
5.99x10-2 
8.87x10'! 
1 .20x10® 
-1.08x10"® 


Table Al: Extracted parameters (from the magnetization data of [2]) for 
the concentration of ATSs and Fe impurities of the BK7 63.530 

g mol~^ ID). The vertical offset represents the residual Larmor diamagnetic 
contribution. 
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BK7 glass (B = 3 T) 



Figure Al; The best fit of the magnetization data [5] for BK7, using Eq. m 
(for the Fe^’*' and Fe^’*' impurities) and Eq. (1191) (for the ATSs). Data from [2]. 
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